arXiv:astro-ph/9706054vl 5 Jun 1997 


The Maximum Optical Depth Towards Bulge Stars Prom 
Axisymmetric Models of the Milky Way 

Konrad Kuijken^ 

Kapteyn Instituut, PO Box 800, 9700 AV Groningen, The Netherlands 


ABSTRACT 

It has been known that recent microlensing results towards the bulge imply mass 
densities that are surprisingly high given dynamical constraints on the Milky Way mass 
distribution. We derive the maximum optical depth towards the bulge that may be 
generated by axisymmetric structures in the Milky Way, and show that observations 
are close to surpassing these limits. This result argues in favor of a bar as a source 
of significantly enhanced microlensing. Several of the bar models in the literature are 
discussed. 
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1. Introduction 


Large surveys for microlensing events in the Galaxy 
have grown into a serious tool for evaluating models 
of the mass distribution of the Milky Way, as is evi¬ 
dent from the recent results of the MAC HO ( Alcock 
t al. 19971), OGLE ([Udalski et al 19941), EROS (|Re¬ 


[lault et al. 1996) and DUO ( Alard, Mao fc Guib- 
rt 1995|) surveys. Microlensing occurs whenever the 


light from a distant star passes sufficiently close to an 
intervening object (the ‘lens’) that its path is signifi¬ 
cantly deflected. The deflection results in a focusing 
of the source star’s rays, which may be detected as an 
apparent brightening of the source with a character¬ 
istic light curve. Typically the brightening lasts for 
a number of days or months before the source and 
lens proper motions spoil the alignment. Both the 
optical depth r, which is the average fraction of the 
sources that is being amplified more than a factor of 
1.34 at any given time, and the average duration of 
the individual events, depend on the distribution and 
masses of leasing objects between the observer and 


the sources (Paczyhski 1986). 


In this paper we concentrate on microlensing to¬ 
wards the bulge of the Milky Way. Both the MAGHO 
and OGLE teams find that the optical depth in this 
direction is rather higher than had been expected ini¬ 
tially: the optical depths reported are in the range 
of 3-6x10“®, roughly twice the expected value based 
on standard Galaxy models. Although initially this 
discrepancy was attributed to our Galaxy’s central 


bar ([Paczyhski et al. 1994 

, Zhao, Spergel & Rich 

1995, Zhao, Rich & Spergel 1996 

), evidence for which 


had been accumulating for some time (e.g. Blitz & 


Spergel 1991, Weinberg 1992; for a review, see Kui 


jken 1996 ), the most recent models for the bar in fact 
do little to enhance the expected optical depth over 
that expected from axisymmetric models (Bissantz et 
al. 1996 , Nikolaev fc Weinberg 1997 ). It therefore 
seems opportune to return to the axisymmetric mass 
models once more. 


In this letter we derive a strong limit on axisym¬ 
metric models: a given optical depth towards the 
bulge comes with a minimum contribution to the 
Galactic circular speed at the sun. It turns out that 
the current values for the optical depth towards the 
bulge are already close to implying a circular speed 
higher than that observed. 

The only result from microlensing theory that is 
relevant for this discussion is the optical depth t to 


a source at distance L, due to a distribution of lenses 
with mass density p{D) at distance D down the line 
of sight: 


47rG 


T = 


P{D) 


D{L - D) 


dD. 


( 1 ) 


T is independent of the mass function f{m) of the 
lenses, though the rate and detectability of microlens¬ 
ing events are not: the survey teams correct for de¬ 
tection (in)efficiency in their analyses. 


2. Can Axisymmetric Mass Components Ex¬ 
plain Bulge Microlensing? 

Axisymmetric density distributions can be viewed 
as constructed from concentric rings. We therefore 
start by considering the relation between the mass 
and the optical depth of a ring. We assume that the 
sources are on the Galactic minor axis, at a distance 
of Ro = 8kpc, and not mixed spatially with the ring; 
these results are thus appro priate for the micro lensing 
of the bulge clump giants ( [Alcock et al. 1997 ). 


2.1. Optical Depth from a Ring 


Consider a ring of uniform surface density with ra¬ 
dius R < Ro and radial width AR R. Its optical 
thickness towards bulge sources near the Galactic mi¬ 
nor axis, at latitude b radians (assumed small), is (cf. 
eq. 0) 

_ 47rG (Ro -R)R 


Ro 


-p{z)AR. 


( 2 ) 


Here p is the density where the line of sight intersects 
the ring, at height z = {Ro — R)b from the plane. 
Now, if the ring is axisymmetric with uniform surface 
density S, then its mass is 


Mh 


2ttRARTi = 


C 'TringRo^ 

2G 



(3) 


(using Ro — R = z/b). The bracketed factor is dimen¬ 
sionless, and depends only on the shape of the vertical 
density profile of the ring, not on the total mass of the 
ring. We will call it the vertical shape factor F^. Its 
denominator tends to zero at small and large z (pro¬ 
vided the ring has finite surface mass density), and 
therefore this geometrical factor has a minimum at 
some intermediate z. Hence equation ^ can be used 
to derive a lower limit on the mass of an axisymmetric 
ring of a given optical depth r. 
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The lowest vertical shape factor for monotonic 
p(\z\) is obtained when the density is constant up 
to a height a > z, and zero beyond. In this case 
E = 2ap, and > 2. More realistic density profiles 
yield values considerably larger, though. Table lists 
the lower limits to for some different forms of p(z). 
In light of those results, we adopt F^ ^ 4.0. 


the lower mass limit of Eg. 3 can be recast as a 


cir cular velocity at the solar circle: approximating 


models for the disk mass distribution and raising their 
mass until their rotation curve reaches 200km s“^ at 
any point. We consider four functional forms for the 
disk surface mass density, whose rotation curves are 
shown in Fig. |^. We will also compare these models 
to solar-neighborhood normalized models, where the 
surface mass density in stars is ~ SSMqpc”^ (Kuijken 
fc Gilmore I989| ). Note that the total surface density 
of matter within I kpc of the plane, 7IM0pc“^ (Kui- 


the circular speed due to a ring by its monopole term jken fc Gilmore 1989| ) is consistent with such maxi- 


(itself a lower limit), we obtain 

.b 


I’oiTc > 


GMring ^ 2_ 

D — o ^^ /ring 

/to Z 


b. ( 4 ) 


Since the radius of the ring is no longer explicit in 
this equation, this limit applies equally to superposi¬ 
tions of rings, i.e., axisymmetric mass distributions. 
Scaling to observed parameLers, we ublaiii - 


^^circ > 




4.0 


( 5 ) 


as the minimum circular speed generated at the solar 
pos ition by axisymmetrically distributed matter re — 


spo nsible for an optical depth of r_o x 10 ^ towards 


bulge sources at Galactic latitude b. A radially ex¬ 
tended mass distribution, in order to lens optimally, 
needs to have F^ close to its minimum at all radii: 
this is achieved when the scale height of the mass is 
proportional to i?o — i?, i.e., tapers linearly to zero at 
the solar radius. 

2.2. Specific Disk Models 

We know more about the Milky Way rotation curve 
than its value at the solar position: we also know that 
the circular speed does not greatly exceed the local 
value at any point in the inner disk of the Galaxy. We 
now investigate how this extra constraint affects the 
maximum disk optical depth, by considering various 


Table I: The minimum vertical shape factor F^ 
Ti/(zp) for different vertical density profile shapes. 
p{z) Attained at z = 


constant to z = a 

2 

a 

exp(-|z|/a) 

2e 

a 

sech(z/a) 

4.7 

1.2a 

sech^(z/a) 

4.5 

0.8a 

exp(—z^/a^) 

4.1 

0.7a 


mal disk models only f or short 2.5kpc) scale length 
models ( Sackett 1997|) . 

The traditional disk model is a radially exponen¬ 
tial disk, E(i?) oc exp(—We consider two verti¬ 
cal density profile shapes, a vertical exponential and 
a sech^(z/a) ‘isothermal’ model. If the scale height 
is small, the rotation curve peaks at a radi us of 2.2 
scale lengths at the value 0.62(GM/h)^/^ ( Freeman 
I970| ). Figure || shows the optical depth from such 


disks towards the Galactic minor axis at a latitude of 
b = 3.5°, for various scale heights and lengths. 

For comparison, we consider also two variations on 
the exponential disk: a less centrally concentrated 
model with S(i7) oc Rexp{—R/h) (e.g., Caldwell & 


Ostriker 1981), and a more concentrated one obtained 


by adding a second exponential disk of a quarter the 
scale length, and one eighth the mass of the main 
component. As may be seen from Figs. ^ and the 
maximum optical depth from such disks is compara¬ 
ble to the exponential disk. 

Finally, we consider a Mestel disk ( jMestel 1963| ), 
which has surface mass density E(i?) = / (27rGi?), 
and whose circular speed is everywhere equal to V. 
Such a model microlenses more efficiently than an ex¬ 
ponential disk since it can satisfy our rotation curve 
constraint at every radius simultaneously. The max¬ 
imum optical depth from a Mestel disk towards b = 
3.5° turns out to be 2.5 x 10“® for a vertically sech^(z) 
mass distribution. 

We summarize the results for the various disk mod¬ 
els in Table |[ We give two numbers for each model: 
the maximal optical depth for each model assuming 
constant scale height, and the (higher) optical depth 
that can be achieved by tapering the disk scale height 
down to zero at the solar position, in proportion to 
(i?o - R)- 
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Fig. 1.— The rotation curves for three types of disk 
models considered, in units of {GM. 



Disk Scale Length (pc) 


Fig. 2.— The optical depth towards sources at the 
Galactic Minor axis from radially exponential disks 
with maximum circular speed 200kms“^, and differ¬ 
ing scale heights and lengths. The solid lines refer to 
vertically exponential disks, the dotted ones to sech^ 
disks. Sun symbols show the optical depths from disks 
normalized to the ‘locally identified’ stellar surface 
density of 35 Mqpc“^ at the solar position, vertically 
exponential with scale height 300pc. 



Fig. 3.— As Figure ^ but with a central hole in the 
disk mass profile. 



Fig. 4.— As Figure but with a more concentrated 
mass profile. 


4 































3. Discussion 


The calculations presented above show that ax- 
isymmetric models of the Galaxy cannot avoid high 
circular velocities at the sun if the optical depth to¬ 
wards the Galactic minor axis is as high as is being 
reported. The current numbers are in fact so high 
as to exceed the maximal disk values derived in Ta¬ 
ble They are also in conflict with the cl aimed de¬ 
tec tion o f the dark halo toward s the LMG ( Alcock et 
al. 1996 , Renault et al. 1996 ), since the small disk 
scale heights 300pc) required for bulge microlens- 
ing would produce virtually no microlensing towards 
the LMG. If indeed the optical depth is too high, then 
non-axisymmetric models, particularly the Galactic 
bar, should be considered. 

Several in-depth studies of microlensing by the cen¬ 
tral bulge/bar have been carried out already. While 
the recent investigations by the Oxford group (Binney 
et al. 1996|; pdissantz et al. 1996 ) and by Nikolaev & 
We .nberg 1997 conclude that the enhancement in op¬ 
tical depth over a standard bulge model ( Kent 199^ ) is 
small, the earlier analysis by Zhao et al. (1995,1996) 
did find a significant enhancement in r. These dif¬ 
ferences illustrate the difficulty in disentangling the 
structure of the inner Galaxy from our view point in 
the disk. 

The Zhao et al. and Oxford studies are based on 
the COBE/DIRBE map of near-infrared emissions, 
the asymmetries of which are used to infer the bar’s 
orientation and shape from perspective effects. Zhao 
et al. ba sed their model on the fitting to the DIRBE 
data by Dwek et al. 199"^ , whose best-fit model is 
quite elongated, with major axis close to the line of 
sight. This model is therefore efficient at microlens¬ 
ing towards the bulge, since it places most of its 


Table 2: Maximum optical depth towards stars on 
the Galactic minor axis at latitude b = 3.5°, for var¬ 
ious disk surface density profiles E(i?). All models 
have been normalized to have a peak circular speed 
of 200kms“^, and have vertical density ocsech^)^). 

The optical depths for models whose scale height is 
radially constant, and tapered oc {Rq — R) are shown. 

Disk T.{R) Tconst X 10^ rtapered X 10^ 


oc exp(—i?/a) 

2.0 

2.6 

oc Rexp{—R/a) 

1.9 

2.4 

oc exp{ — R/a) + 2 exp( —i?/4a) 

2.1 

2.7 

oc R~^ 

2.5 

3.3 


mass density between us and the bulge. The Ox¬ 
ford group, on the other hand, use a Lucy deprojec¬ 
tion of the emission to infer the bar properties. They 
find a much shorter bar, which consequently gener¬ 
ates less microlensing. Nikolaev & Weinberg derive 
their bar parameters from star counts of IRAS AGB 
sources. Their bar is considerably longer than the 
GOBE-derived ones, but it is quite broad, providing 
only a modest enhancement in lensing over axisym- 
metric models. 

It is difficult to prefer one of these models over an¬ 
other. The Zhao et al. model is a simple parametric 
fit to the GOBE data, whereas the Oxford model is 
completely non-parametric. In fact, the Oxford inver¬ 
sion is puzzlingly robust: different starting conditions 
result in very similar solutions, in spite of the fact that 
a three-dimensional emissivity model is reconstructed 
out of a two-dimensional projected surface brightness 
map! One assumption made is, however, crucial: the 
models are taken to be 8-fold symmetric, i.e., the bar 
is modeled as reflection symmetric in three orthogonal 
planes (it is this symmetry that allows the observed 
left-right differences in the map to be interpreted as 
perspective effects). In reality, even the most regu¬ 
lar bars observed do not exhibit this symmetry: bars 
are generally somewhat twisted, and often connect 
onto spiral arms joining at the ends. As Binney et 
al. themselves recognize, even a 4-fold symmetry (re¬ 
flection in Galactic plane and through the Galactic 
center) is insufficient to allow a stable inversion. 

The Weinberg bar may be stronger than it appears 
in the analysis of IRAS sources, since a relatively 
small number of harmonics was used in modeling the 
star counts (m = 0,2,4 azimuthal harmonics). It is 
therefore possible that the true density enhancement 
towards the bulge is higher than modeled. This bar 
is also clearly larger than the GOBE models, suggest¬ 
ing that the Galaxy may well be a bar-within-a-bar 
galaxy, a type of system which is relatively common. 
In that case, the GOBE models may need to take ac¬ 
count of this extra foreground component. 

4. Conclusion 

We have shown that if the mass responsible for the 
microlensing towards the Galactic bulge is distributed 
axisymmetrically, the high optical depths r being re¬ 
ported may be translated into simple but significant 
lower limits on the circular speed at the solar position 
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due to lensing matter: 

v1 > 2.0c^r6 (6) 

where b is the galactic latitude (in radians) of the 
sources. This limit evaluates to Vc > 210kms“^ for 
an optical depth of 4 x 10“® at 6 = 3.5°, very close to 
the actual rotation speed of 200±20kms“^ 

19c(2|). 

This limit is approached when the lensing mat¬ 
ter has a scale height that is well matched to the 
height of the line of sight towards the stars being 
monitored, otherwise even higher circular speeds are 
implied. More restricted, but realistic, models for the 
mass distribution fail to generate optical depths above 
2.5x10“® without requiring circular speeds in excess 
of 200km s“^. 

This result comes close to ruling out purely ax- 
isymmetric models for the mass distribution in the 
inner Galaxy. Those models left require a very flat¬ 
tened distribution of the Galactic mass, of thickness 
^ 300pc, contrary to the conclusions drawn from the 
microlensing experiments towards the Large Magel¬ 
lanic Gloud, which favor a round dark matter distri¬ 
bution. A significant enhancement in r from a central 
bar appears to be a more plausible explanation. 

It is a pleasure to thank Penny Sackett for discus¬ 
sions and a careful reading of the manuscript. 
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